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Abstract The structure function is a useful quantity to characterize wavefront distortions. 
We derive expressions for the structure functions of the averaged wavefront phase and slopes. 
The expressions are valid within the inertial range of atmospheric turbulence, and are meant 
to serve as engineering formulae when reconstructing profiles of the atmospheric turbulence, 
specifically in the context of atmospheric profiling instruments (e.g. SLODAR and S-DIMM+) 
and multi-conjugate adaptive optical systems. 


1 Structure function 

Kolmogorov’s theory of turbulence provides the structure function of phases in a wavefront. It 
equals the expected mean squared difference of the phases, / and /', at two points, (x, y ) and 
(. x',y '), separated by distance s: 

/ \ 5 / 3 

sF(s) = {(f-m = k\j- j (i) 

with k = 6.88. ro is the Fried parameter: in a circular region with a diameter equal to the 
Fried parameter, ro, the phase variance is roughly equal to 1 (see for example Roddier ED)- The 
symbol {•} stands for the statistical expectation. 

SF(s ) relates to the point function, f(x, y), and to a coherent wavefront. In actual adaptive- 
optics (AO) conditions one deals with averaged phases or phase-slopes, since measurements 
utilize sensors of finite size. Accordingly, the structure functions need to be adjusted to apply 
to the averaged phases or slopes. 

Expressions for the adjusted structure functions have been derived by a number of authors, 
see Tokovinin [2] and references therein. These expressions were derived for use with the Dif¬ 
ferential Motion Monitor (DIMM) [3] and are valid for separations larger than the averaging 
diameter, s > d. 

In order to reconstruct profiles of the atmospheric turbulence, it is useful to extend the ex¬ 
pressions to separations that are smaller than the averaging diameter, see for example Scharmer 
& van Werkhoven [T] and Kellerer et al. [5]. Here we derive expressions valid within the inertial 
range of atmospheric turbulence, i.e. the range of separations where Eq.|T| correctly approxi¬ 
mates the structure function of the phase. 
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2 Structure function of the degraded phase 


Let g be the phase averaged over a (circular or square) region A centered at (x,y): 

g(x,y) —< f(x',y') > A (2) 

The notations adopted in the following text are explained in AppendixlAl 
As with the undistorted phase, /, the variance of g cannot be given without specifying a 
reference region. This is not the case for the structure function. 

Let A be a square of size dx d or a circle of diameter d. For large distances s the structure 
function, SF(s;d), of g converges to the structure function, SF(s), of /. How the functions 


differ at smaller distances s needs to be investigated: 

SF(s;d) = {(g(x,y) - g(x + s,y)) 2 } (3) 

= {(< f(x,y) >a - < f{x + s,y) > A ) 2 } (4) 

= 2{< (f(x, y)>\-< /(A y) >a < f(x + s, y) >^} (5) 

= 2 {< f(x, y)f(x ', y ') - f{x, y) f(x' + s, y) >, 4 } (6) 

= {< (/Or, y) - fix' + a, y')f - (f(x, y) - f(x f , y')) 2 ) > A } (7) 

= k < (u ,5/3 - u 5/3 ) >^4 /r 0 5/3 ( 8 ) 


with: u = ((x — x') 2 + (y — y') 2 ) 0 ' 5 and v! = ((x — x' + s) 2 + (y — y') 2 ) 0 ' 5 . In standard notation 
this result reads: 


SF(s ; d) = k [ [ (u ,5/3 - u 5/3 ) dx dy dx’ dy'/S 2 (9) 

J A J A 

S: area of A ( d 2 for a square, 7rd 2 /4 for a circle). 

To put this into words: Let A' be the averaging region A shifted by distance s in the 
direction x. The degraded structure function is then equal to the average structure function 
between a point in A and a point in A', minus the average structure function between two 
points within A. For s » r the degraded structure function, SF{s]d), converges to SF(s). 

Since both SF(s;d ) and the unmodified structure function SF(s ) = k(s/ro ) 5 ^ 3 increase 
steeply with s the comparison of the structure function can best be made in terms of the ratio 
of the two functions. This ratio can be termed the reduction factor, 


RF(s] d) 


SF(s-,d ) 
k (s/r 0 ) 5 / 3 


( 10 ) 


It specifies - for two points separated by distance s - the reduction of the mean squared phase 
difference due to the phase averaging over diameter d. 


Approximate formulae The numerical evaluation is shown on Fig. CD it is obtained from 
10 6 point pairs chosen randomly within a circle of diameter 1. The following analytical function 
approximates the result with less than 0.5% deviation over the range s/d = 0.01 — 100 (see 
right panel of Fig.[l]): 

RF(s ; d) ~ (1 + 1.14 (s/d)- 5 ' 5/3 )- 1/5 ' 5 
/ \ 5 / 3 

SF(s ; d) ~ k (! + 1-14 (s/d)~ 5 - 5/3 )- 1/5 - 5 
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Figure 1: Left panel: Reduction function obtained from 10 6 random point-pairs. The standard 
deviation is too small to be identified. Dots: Numerical evaluation, line: analytical approx¬ 
imation in terms of Eq.[[2j Right panel: Deviation of the analytical approximation to the 
numerical values. 


3 Structure function of the wavefront slope 

Line-averaged slope The simplest approximation to the slope is the phase difference be¬ 
tween two points divided by their distance, d: 

s x(x, y) = 7 ^ (/ (x + d,y) — f(x, y))/d (13) 

s y (x, 3/) = ^: (/(^ V + d )~ f( x » y))/ d ( 14 ) 

Mean slope A related notion is the slope averaged not over a line element, d, but over a 
reference region, A. It equals the derivative of the degraded phase g(x,y ) (see Eq.[2]): 

s x =< d f(x', y')/dx > A = d g(x, y)/dx (15) 

s v =< d /(^> y')/&y >A= d g(x, y)/dy (16) 

Least-square slope A further concept is the slope of the least-square linear approximation 
to the wavefront over the reference region. It is particularly relevant because it is a close 
approximation of the quantity measured by a SH-sensor [7]. 

To avoid an excess of symbols, the same letters, s x , s y , are used for the three different 
choices; the type of average needs to be recognized from the context. 


3 





3.1 Average gradient over a line segment 

The simplest approximation to the slope is the phase difference between two points a distance 
d apart: 


s x{x, y) = (f(x + d,y) — f(x, y))/d (17) 

8 y (?, y) = ^ (/Od y + d )~ f ( x > y))/ d ( 18 ) 

In the equation for the structure function of s x only the variable x appears, which makes 

it convenient to use the simpler notation f x for f(x,y), etc. In this simplified notation the 
structure function of the slope, s x , in the direction of the distance x, reads: 

SF x (s ; d) = {(s x (s, 0) - s x (0, 0)) 2 } (19) 

= ( tx ) {(f s+d ~ ~ f d + ( 2 0 ) 

= (T) ‘ -/„) 2 } + {</,- -{(u-f,f} (2i) 

- {(/o - /,«) 2 } + {(/o - /,) 2 } - {(/,« - /,) 2 })/<i 2 (22) 

The classical DIMM formulae: If s » d, this simplifies to: 

SF x (s-, d) = ^ (2 SF(d)/d 2 - SF"(s )) (23) 

= 2k (^j r^ 5/3 (d _1/3 - 5/9 s" 1/3 ) (24) 

= 0.35 A 2 r“ 5/3 d _1/3 (1-5/9 (s/d)~ 1/3 ) (25) 

Similarly, for the transverse slopes: 

SF y (s] d) = 2 (SF(d) + SF(s) - SF((s 2 + d 2 )°' 5 ))/d 2 (26) 

With 

SF{(s 2 + d 2 ) 0 ' 5 ) = k(s 2 + d 2 ) 5/6 r" 5/3 (27) 

~ k (s 5/3 + 5/6 s- 1/3 d 2 ) r" 5/3 = SF(s) + 5/6 k s" 1/3 d 2 r" 5/3 (28) 

one obtains the equivalent of Eq.l25l 

SF y (s ; d) = 0.35 A 2 r" 5/3 d~ 1/3 (1 - 5/6 (s/d)~ 1/3 ) (29) 

The structure functions of the x— and y— slopes - as used for DIMM - are represented as blue 
dots on Figs. [2] and [3j 
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General relations: Eqs.Hni and [2H] are limited to separations larger than the reference diam¬ 
eter: s » d. They are the classical relations indicated by Sarazin and Roddier for the analysis 
of DIMM measurements. More general relations, valid for any value s, are: 

SF x (s ; d) = k ^ r" 5/3 (2 d 5/3 - \s - d| 5/3 + 2s 5/3 - (s + df/ 3 )/d 2 (30) 

= k (J-J r^d - 1 ' 3 (2 - |1 - s/d| 5/3 + 2 (s/d) 5 / 3 - (1 + s/df /3 ) (31) 

SF y {s- d) = 2k (J-'j 2 r~ 5/3 (d 5 / 3 + s 5 / 3 - (s 2 + d 2 f/«)/d 2 (32) 

= 2k f A^ r~ 5/3 d - 1 / 3 (1 + ( s/df 13 - (1 + (s/d) 2 ) 5/6 ) (33) 

These two structure functions are indicated as blue lines on Figs.[2]and[3] 


3.2 Average gradient over a reference region: G-tilt 

A related notion - which applies to a defocussed wavefront image - is the slope averaged not 
over a line element, s, but over a reference region, A. It equals the partial derivative of g(x,y ) 
(see Eq.|2]): 

=< d f(x', y')/dx > A = dg(x, y)/dx (34) 

s y =< df(x', y')/dy > A = d g(x, y)/dy (35) 

The relevant structure functions can be derived from the structure function of the averaged 
phase: 

SF x (s ; d) = liniA^o (2 SF( A; d)/ A 2 - SF"{s ; d)) (36) 

SF y {s- d) = 2 li mA ^ 0 (SF{ A; d) + SF(s ; d) - SF({s 2 + A 2 ) 0 ' 5 ; d)) /A 2 (37) 


Approximate formulae We use the approximation for the structure function of the de¬ 
graded function (see Eq. [T2]) to derive the expressions for the structure function of the gradient 
over a circle: 

SF X ( S] d) = f A V k d- 1/3 r~ 5/3 (2 a- 1 ^ - ± u^ 3 (10 + a (7 - b) u~ b ' 3 w~ l + a 2 (b + 1) u 2b ' 3 w~ 2 ) 

\2ir J V 9 

(38) 

SF y (s;d) = 2 (F\ kd,- 1/3 rF 3 ( a" 1 '*- ) o" 1,J (5 + uu-*' 3 ui"i)j (39) 

where a = 1.14, b = 5.5, u = s/d and w = 1 + au ~ b A This rewrites as, 

SF x (s ; d) = 2 k d" 1/3 r~ 5/3 (l.95 - i u ~ 1/3 (10 + 1.71 u ~ 5 ' 5 / 3 w~ l + 8.45 u ~ 11/3 «T 2 )) 

(40) 

/ ^ ) kd-^r- 5 ' 3 (o.dS-^u-^w-^ib + l.Uu-^w- 1 )') (41) 


SF y (s-d) = 2 


X 

2tt 
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These structure functions are represented as black dots on Figs. [2] and El The functions are in 
excellent agreement with the expressions given in Tokovinin [2] for s > d: see his Eq. 7. Note 
that Tokovinin’s equations are based on calculations by Conan et al. [6], 


3.3 Least square fit to the wavefront: Z-tilt 

Another concept is the slope of the least-square linear approximation to the phase front over 
the reference region. This variable will be considered, because it equals the tilt measured by a 
Shack Hartmann (SH)-sensor: 

Let (0, 0) be the center of the circular reference domain, C, of diameter d. To compute the 
tilts in x— or y —direction over C, one minimizes the mean squared deviation, S, between the 


particular phase pattern, and a linear approximation: 

S =< (a 0 + aix + a 2 y - f(x,y)) 2 > c (42) 

< .. > stands for the normalized integral, i.e. the mean value, over C. 

dS/da 0 = 0 = a 0 + ai < x >c +a 2 < y >c ~ < f(x, y) >c (43) 

dS/dai = 0 = a 0 < x >c +a± < x 2 >c +a 2 < xy >c — < x f(x, y) >c (44) 

dS/da 2 = 0 = a 2 <y > c +eq < xy > c +a 2 < y 2 > c - < y f(x,y) > c (45) 

Due to the symmetry of the reference domain the moments < x >, < xy > and < y > vanish. 

< x 2 >c=< y 2 >c= qd 2 , with q = 1/16 for a circular and q = 1/12 for a square domain. Thus: 

a 0 =< f >c (46) 

ai =< x f(x , y) > c /(q d 2 ) (47) 

a 2 =< y f(x, y) > c /(q d 2 ) (48) 

If s is parallel to the direction of the slope (ai), the structure function equals: 

AF x ( S ;d) = {(a 1 (0)-a 1 ( S )) 2 } (49) 

= {(< x f(x, y) > c — < x f(x + s, y) >c) 2 }/(q 2 d 4 ) (50) 

= 2{< (x f(x, y) > 2 C - < x f(x, y)> c <x f(x + s, y) > c }/(q 2 d 4 ) (51) 

= 2{< x x' f(x, y) f(x', y') - x x f(x, y) f(x' + s, y') >c}/{q 2 d 4 ) (52) 

= {<xx' (f(x, y ) - f(x' + s, y')) 2 -xx' (f(x, y ) - f(x', y')) 2 ) > c }/{q 2 d 4 ) (53) 

= k r o 5/3 < —xx'(u' 5 / 3 — u 5 / 3 ) >c /(2 q 2 d 4 ) (54) 

with: u — {(x — x') 2 + (y — y') 2 ) 0 ' 5 and: u' = ((x — x' + s) 2 + (y — y') 2 ) 0 ' 5 . 

Similarly, if s is perpendicular to the direction of the slope (a 2 ), the structure function 

equals: 

SF y ( S] d) = ^ k r~ 5/3 < -yy\u^ 3 - u 5 / 3 ) > c /(2 q 2 d 4 ) (55) 
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Figure 2: The structure function of the x-slope. Solid black line: Least square slope over a 
circle of diameter d (Z-tilt). Black dots: Average slope over a circle of diameter d (G-tilt). Blue 
line: Average slope over a line segment of length d. Blue dots: Approximation used for the 
DIMM. Left panel: log-log scale, right panel: log-linear scale. 


Approximate formulae The numerical evaluation is shown on Fig.[H it is obtained from a 
two-fold integration over a 60 x 60 grid. The following equations approximate Eqs.lMl and Pol 
with less than 3% deviation for 0.01 < u = s/d < 100: 

SF x (s] d) = k d~ 1/3 r~ 5/3 (2.06 - 1.55 (1 + 10.5 u 2 )~ 1/6 - 0.51 (1 + 10 u 3 - 3 )~ 2/3 - 3 ) (56) 

SF y (s] d) — kd~ 1/3 r~ 5/3 ( 2 .O 6 -1.72 (1 + 1.5 w 2 )~ 1/6 - 0.34(1 + 1.4 m 2 )' 1 ) (57) 

The numerical evaluation are represented as black lines on Figs.[2]and[3] The functions are in 
excellent agreement with the expressions derived by Tokovinin [2] (and based on calculations 
by Sasiela [8]) for s > d: see Eq. 8 in Tokovinin [2]. 


4 Conclusion 

We have derived structure functions for the averaged phases and slopes from a 5/3 power law 
for the structure function of the phase. The expressions are summarized in Tabled and are 
valid within the inertial range. For larger separations, the phase structure function is given by 
(see Eq. 3.22 in Conan [9]): 


SF(s-L 0 ) 


2 1 /er(ll/6) 





2 v 6 
UV6) 



(58) 
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Figure 3: The structure function of the y-slope. Solid black line: Least square slope over a 
circle of diameter d (Z-tilt). Black dots: Average slope over a circle of diameter d (G-tilt). Blue 
line: Average slope over a line segment of length d. Blue dots: Approximation used for the 
DIMM. Left panel: log-log scale, right panel: log-linear scale. 






s/d s/d 

Figure 4: Left panel: Structure function of tip (black) and tilt (blue) over a circle of diameter d. 
The dots show the numerical evaluation, the lines are analytical approximations (Eqs. [56H57|) . 
The standard deviation is too small to to be identified. Right panel: Deviation of the analytical 
approximation to the numerical values. 
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Quantity 

Structure function 

Phase 

k (s/r 0 ) 5/3 (1 + 1.14 (s/d) -5 ’ 5 / 3 ) -1 / 5 ' 5 

x-slope 

DIMM approximation 

2(1 - 5/9 u -1 / 3 ) 


Average gradient over a 
line segment of length d 

2 — |1 — ri| 5 / 3 + 2u 5 / 3 - (1 + u ) 5/3 


Average gradient over a 
circle of diameter d (G-tilt) 

1.95-1 /9 u~ 1 / 3 n; -1 / 5 ’ 5 (10 + 1.71 ri -5 ’ 5 / 3 w -1 + 8.45 u -11 / 3 w~ 2 ) 


Least square slope over a 
circle of diameter d (Z-tilt) 

2.06 - 1.55 (1 + 10.5 u 2 ) -1 / 6 - 0.51 (1 + 10 «3-3)-2/3.3 

y-slope 

DIMM approximation 

2(1 - 5/6 u -1 / 3 ) 


Average gradient over a 
line segment of length d 

2(1 + u 5/3 - (1 + it 2 ) 5 / 6 ) 


Average gradient over a 
circle of diameter d (G-tilt) 

2(0.98 - 1/6u -1 / 3 u> -1 / 5 ’ 5 (5 + 1.14tt -5 ’ 5 / 3 u -1 )) 


Least square slope over a 
circle of diameter d (Z-tilt) 

2.06 - 1.72 (1 + 1.5 u 2 ) -1 / 6 - 0.34 (1 + 1.4 it 2 ) -1 


Table 1: Structure functions for the degraded phase and slopes, u = s/d and w = 1 + 
1.14 u -5 ’ 5 / 3 . The formulae for the x- and ^/-slopes need to be multiplied by an additional factor 
(A/27 t) 2 k d -1 / 3 Tq 5 ^ 3 . 

T: Gamma function, K 5 / 3 : Bessel function of the third kind and of order 5/3. A good 
approximation of this expression, is: 

SF(s ; L 0 ) ~ k ^1200 + 60 

with k = 6.88. Fig. 0 compares the exact expression (Eq. 158]) with the analytical approximation 
(Eq. l59j) and shows that the the 5/3 regime breaks off for separations s > Lq/100. The present 
approximations are thus meant to serve as useful engineering formulae when dealing with 
small separations, e.g. when measuring profiles of the atmospheric turbulence with site-testing 
telescopes such as SLODAR j IT) and SDIMM+ [T]. 

A Shorthand notation for expectation values of phase 
integrals 

In the text certain averages over the reference domain, A, are considered, that are integrals of 
the phase, over A, or are related quantities, such as the product of the phases, f(x,y ) 

and f(x',y'), of all point pairs within the region. 

To make the equations more transparent, a shorthand notation is used for the integrals. For 
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Figure 5: Left panel: Structure function of the phase with finite outer scale, L 0 . The dots 
represent the exact values as calculated by Conan |9j (Eq.[5H]), the black line shows the analyt¬ 
ical approximation fEq. l^Ul) . The blue line indicates the structure function with infinite outer 
scale (Eq.[T]). Right panel: Relative difference between the exact values and the analytical 
approximation. 


example: 

< f(x,y) f(x',y') > A = I [ f(x,y) f (x , y 1 ) dx dy dx dy'/S 2 (60) 

J A J A 

where S is the surface of A, and the integration runs over all point pairs (x , y ) and (x\ y') 
in A. 

If the integration runs only over a function, such as f(x,y), of one point, it can, of course, 
be written as a simple integral. But where - in combination with other terms - it is convenient, 
the double integral can nevertheless be retained, i.e. the shorthand notation can remain the 
same: 

<f(x,y)>A= [ [ f (x,y) dxdydx' dy'/S 2 = [ f(x,y)dxdy/S (61) 

J A J A J A 

Since each wave front is given only up to a constant term, expectation values, such as 
{f{x,y) 2 } or {f(x,y)f(x',y')}, are undefined. In the equations they appear, therefore, only 
in combinations were the undefined terms combine to a sum of differences. In particular 
{f(xi,yi) 2 } — {f(x 2 , 1 / 2 ) 2 } = 0 is used to express the expectation of a sum of phase prod¬ 
ucts in terms of squared differences, i.e. in terms of the structure function: 

2{f(x 1 ,y 1 ) f(x 2 ,y 2 ) ~ f(x 3 ,y 3 ) f(x Al y A )} = {(f{x 3l y 3 ) - f(x A ,y A )) 2 } - {(/(xi,j/i) - f(x 2 , y 2 )) 2 } 

(62) 

= SF(s 34 ) - SF(s 12 ) (63) 

where S 12 and S 34 are the distances between (x\,yi), (x 2 ,y 2 ) and (x 3 ,y 3 ), (x 4 ,y 4 ), respec¬ 
tively. 
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The adopted notation can be used to show the well known fact that the variance, cr 2 , in A, 
equals half the mean squared phase difference between two points (x, y ) and ( xy') in A: 


= {< f(x, V ) 2 >a - < f(x, y) >\} (64) 

= {< f(x, y) 2 >a < f(x, y) f(x', y') > A } (65) 

= {<(f(x,y)-f(x',y')) 2 > A }/ 2 (66) 

=< SF(s) > A /2 (67) 

= k/2 < s 5/3 > A /r 0 5/3 (68) 


with: s = ((x — x') 2 + (y — y') 2 ) 0 5 . For a circle of diameter, d: 


<5 5 / 3 >=0.3 d 5 / 3 


aj = 1.02 



(69) 

(70) 
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